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We apply the fermion renormalization group method [1], implemented numerically in Ref.[2], to a 
two-band model of FeAs-based materials. At half filling we find the (tt, 0) or (0, tt) spin density wave 
order and a sub-dominant superconducting pairing tendency. Due to a topological reason, the spin 
density wave gap has nodes on the fermi surfaces. Away from half filling we find an unconventional 
s-wave and a sub-dominant (1^2 _y2 pairing instability. The former has s symmetry around the hole 
fermi surface but exhibits s + d^2_y2 symmetry around the electron pockets where the 90 degree 
rotation is broken. The pairing mechanism is inter-pocket pair hopping. Interestingly, the same 
interaction also drives the antiferromagnetism. 



Recently there is a flurry of interest in the the com- 
pound MFeAsOi_ajFa;[3]. It has been shown that by 
varying the rare earth elements (M)[4], these materials 
can be made superconducting with Tc as high as 55K[5]. 
This has stimulated a flurry of interest in these materials. 
At the present time, preliminary experimental results 
have indicated that these materials exhibit semimetallic 
antiferromagnetism at stoichiometry (i.e., x = 0). Upon 
substituting O with F the antiferromagnetism diminishes 
while superconductivity appears [6]. 

Structurally MFeAsOi_a^Fa^ can be viewed as As-Fe- 
As tri-layers separated by MOi-x^x spacers. Early 
electronic structure calculations [7] and angle-integrated 
photoemission[10] suggest that the carriers at the fermi 
energy are essentially Fe in character. Thus, like many 
others, we will focus on a As-Fe-As tri-layer in the fol- 
lowing discussions. To envision these tri-layers, imagine 
a square lattice of Fe. The As sit either above or below 
the center of the square plaquettes, and these two types 
of plaquette form a checkerboard. Each structural unit 
cell contains two Fe (dashed rectangle in Fig. (l)(a), and 
the basis vectors are X = f + ^ and Y = x — y where x 
and y are the basis vectors of the Fe square lattice (see 
Fig. (l)(a)). In the following the we shall refer to the 
reciprocal unit cell associated with the basis vector x^y 
and X, Y as the "unfolded" and "folded" Brill! ouin zone 
(BZ) respectively. 

Aside from the semi-metalicity, there are other impor- 
tant differences between the antiferromagnetism in this 
material and that of the cuprates. First, in the unfolded 
BZ the magnetic ordering wavector is either (0, tt) or 
(7r,0)[ll] instead of the usual (tt, tt). Second, the or- 
dering moment is quite small (~ 0.25/iB)[ll] compared 
to that of the cuprates. For the superconducting state, 
penetration depth [12], Hc2[13], 1/Ti of nuclear spins [14], 
/iSR[15], and point tunneling measurements [16] all in- 
dicate the presence of line nodes in the superconduct- 
ing gap. However, more experiments, in particular those 
done on single crystals, will be necessary to check the 



above conclusions. 

On the theoretical side, several electronic structure cal- 
culations suggest the presence of hole and electron pock- 
ets at fermi energy. It is suggested that the near nesting 
of these pockets is responsible for the antiferromagnetism 
at X = 0[7-9]. Estimate of the strength of local coulomb 
interaction suggests that this system is on the border be- 
tween strong and weak correlation [17]. Although LDA 
type calculations suggest all five orbitals of Fe contribute 
to states at the fermi energy [7], this was simplified re- 
cently where only two out of the five bands are kept while 
preserving the electron and hole pockets [18]. Based on 
different degree of simplification of the electronic struc- 
ture, and different approximate treatment of the elec- 
tronic correlation, a number of groups have studied the 
superconducting pairing instability of these material[19]. 
In addition there are a couple of recent attempts to use 
the point group symmetry to narrow down the pairing 
symmetries [20]. The wide spread in pairing symmetry 
concluded from these studies call for a more unbiased 
assessment of the pairing instability. 

In this paper we study a two band model with 
Hubbard-like and Hund interactions. We study the possi- 
bility of electronically induced pairing by performing one- 
loop renomalization group (NRG) calculation[l]. This 
numerical version of this method (NRG) was applied to 
the cuprates by Honerkamp et al[2]. It was shown that 
in the framework of one-band Hubbard model, interac- 
tions which promote dx2_y2 pairing and (tt, tt) antifer- 
romagnetic order are generated at low energies. We be- 
lieve that owing to the weaker correlation, this method is 
much better suited for MFeAsOi-^^F^. The NRG results 
for five band model [7] will be the subject of upcoming 
publication [21]. 

Model Hamiltonian - As Ref. [18] we simplify the elec- 
tronic structure of the As-Fe-As tri-layter by keeping only 
two Fe orbitals: Sdxz and 3dyz. The Arsenic are viewed 
as merely mediating hopping between these orbitals. Due 
to the relative orientation of the Fe and As it is more con- 
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venient to use as basis 3dxz and 3dy^ where x^y and 
X, y are shown in Fig. (l)(a). Our tight binding model 
include nearest-neighbor and next-nearest-neighbor hop- 
pings. From symmetry considerations there are four in- 
dependent hopping parameters ti^t[^t2^t2 as shown in 
Fig. (l)(a). Among them t[ is due to the direct overlap 
of two neighboring Fe orbit als. All the rest three param- 
eters describe hopping mediated by As. In the following 
we shall label the Fe orbitals 3dxz and Sdyz as a = 1, 2. 
The tight-binding Hamiltonian in the unfolded BZ reads 
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^0 = E E ^LsKti^Kks = E E ^n(k)^ik.^„k« 

k,s a, 6=1 k,s n=l 

K{k) = a{k)I + 5,(k)r, + 5,(k)r„ (1) 

where c^ks annihilates a spin s electron in orbital a and 
momentum k, and V^nks is the band annihilation opera- 
tor. The I^Tx^Ty in Eq. (1) are the 2x2 identity ma- 
trix and Pauli matrices respectively. They act on the 
orbital {cixz^dyz) 

space, and Q^(k) — /i + 2t^ cos ky + 
2 cos kx[t[ + (t2 + ^2) cos ky], bxCk) = 2ti(cos kx — cos ky), 
^^(k) = — 2 {t2 — 1 2) sin kx sin ky , and ei,2(k) = a(k) ^ 
^6a;(k)^ + bz(ky. We have checked that when acted 
upon by the element {g) of the point group {C4v) the 
band operator V^aks Vgi^agiis where rjg = ±1. After 
turning on a proper chemical potential /i we get two hole 
pockets around k = (0,0) and (7r,7r) and two electron 
pockets at (0,7r) and (7r,0). In the rest of the paper we 
shall use the following values for the hopping parameters 
h = 0.38el/,t2 = 0.57eV,t'2 = 0. The value of t[ criti- 
cally determines the superconducting gap function, and 
will be discussed later. 

Now we consider local interactions including intra- 
orbital and inter-orbital Coulomb interaction Ui, and /72, 
Hund's coupling Jh and the pair hopping term. When 
summed together they give 

2 

^int = ^{Ul ^ ni,a,T^i,a,i + ^2^i,l^i,2 
i a=l 

^JH['^4ls42s'Cils'Ci2s + (4lT4uC^2!Ci2T + h.C.)]}. 

Here i labels the unit cell, 5, s' =t, i, and rii^a = ^i,a,T + 
ni^a,i is the number operator associated with orbital a. 
The total Hamiltonian H = Hq + H^^t is the starting 
point of our study. 

We have checked that the bare H has no superconduct- 
ing instability for realistic interaction parameters. Thus 
we decide to perform a numerical renormalization group 
(NRG) calculation where high energy electronic excita- 
tions are recursively integrated out. The hope is that 
the effective interaction generated at low energy would 
show the sign of superconducting instability, in the fol- 
lowing we present the result of such a calculation. Tech- 
nical details of the NRG can be found in Ref. [2]. In 




FIG. 1: (a) The in-plane projection of dxz and dyz orbitals 
and the four independent hopping parameters, (b) b as a 
function of k. The red curves are the fermi surfaces. 



brief, we divide the first BZ into N patches, and at each 
renormalization iteration we sum over the five one-loop 
Feynman diagrams labeled (cl)-(c3) in Fig. (2). The es- 
sential complication in the present study is the presence 
of two different bands and four disjoint pieces of fermi 
surface. It turns out that it is easier to work with the 
folded BZ. This is because folding doubles the number 
of bands so that each band has only one fermi surface 
(Fig. (2) (a)) and the BZ patching scheme of Ref. [2] can 
be directly applied. Because of the existence of four band! 
s, we need to keep track of the band indices as well as 
the momenta in the interaction vertex. Thus we need to 
compute 256 x different interaction vertices at each 
renormalization step in contrast to in the single band 
case [2]. Although our RG calculation is performed with 
the folded BZ we shall present our results using the un- 
folded BZ for simplicity. 

The antiferromagnetism and superconducting pairing 
tendency at half filling. At stoichiometry (x=0) the two 
bands derived from the dxz and dyz orbitals accom- 
modate 2 electrons per unit cell. For t[ = the elec- 
tron and hole pockets are perfectly nested by (0, tt), (tt, 0) 
in the unfolded BZ. In a theory like ours, the antiferro- 
magnetism is due to the the above nesting. For t[ 7^ 
the nesting is imperfect. Hence one might expect a full 
spin density wave (SDW) gap for t[ — and a partially 
gapped fermi surface for t'^ 7^ 0. However, as we will 
show in the following, even for the former case there are 
nodes in the SDW gap, and the reason is topological. 

First, we say a few words about how the results are 
obtained. In the calculation we compute the renor- 
malized interaction vertex function y(ki, a; k2, h\ ka, c; d) 
for two incoming electrons with opposite spins. Here 
a,b,c,d = 1,..,4 are the band indices, and ki,..,k3 
are momenta on the Fermi surfaces. From this ver- 
tex function we can extract the renormalized interac- 
tion in the Cooper or SDW channels as follows. For 
the singlet /triplet Copper channel Vf^f a;p,b) = 
V(k, a; — k, a; p, 6; 6) ± y(— k, a; k, a; p, 6; 6); for the 
SDW channel V^^^(k, a, t/; p, c, 6) = -2V(k,a;p + 
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FIG. 2: (a) The division of the first BZ into N (=20) patches. 
The color curves denote the fermi surface(ti = —0.05). (b) 
The strongest inter-pocket pair tunnehng(for t[ = 0). As in- 
dicated by the red and blue arrows, a ±q pair on one pocket 
fermi surface are scattered to a ±k pair on another pocket. 
Interestingly, the same scattering also drives the antiferro- 
magnetism. This can be seen by noting that the momen- 
tum transfer between the incoming spin up (red) and outgoing 
spin down (blue) electrons is exactly the nesting wavevector 
(tt, 0) or (0,7r) as indicated by the dashed line. (cl)-(c3) The 
summed Feynman diagrams. Here the arrowed solid denote 
Green's function and the dashed lines represent the renormal- 
ized interactions. 



Q, b] p, c; (i), where Q is the ordering wavevecctor, band 
indices 6 = c ± 1 mod 4 and d = a ± 1 mod 4. Since 
k and p only takes N discrete values, we can treat Vff 
and V^^^ as matrices. The few lowest eigenvalues as a 
function of the RG evolution (panel (a)) and their final 
associated eigenvectors (panel (b)) of these matrices are 
plotted in Fig. (3). 

In panel (a) of Fig. (3) we show the N=20 RG evolu- 
tion of the scattering amplitudes in the SDW and two 
different types of superconducting pairing channels as a 
function of aln(Ao/A). Here a = — l/ln(0.97), Aq is the 
initial and A is the running energy cutoff. The bare inter- 
action parameters used are Ui = 4.0, U2 = 2.5, Jh = 0.7 
eV. Clearly as Aq/ A increases the interaction that drives 
SDW (black dots) grows in magnitude the fastest. The 
form factor fsDwO^) associated with the SDW order 

{AsDW,a,b = Ek/5DH-,a,5(k)V^5k+(7r,7r)T^aki' b = a±l 

mod 4) is shown in panel (bl) of this figure. Interest- 
ingly despite the perfect nesting of the fermi surfaces (see 
Fig.2(b)) there are two nodes ! 

To understand the origin of these nodes we observe 
that the top of the hole pockets (situated at k = (0,0) 
and (tt, tt)) is doubly degenerate. According to Berry [22] 
the band eigenfunctions must exhibit non-trivial phase 
as k moves around these degenerate k points. Omit- 
ting the identity term, the i^(k) in Eq. (1), i^(k) = 
bz(k)rz + bx(k)Tx^ is that of a spin 1/2 in a k-dependent 



magnetic field. A plot of b as a function of k for 
^1 = ^2 = is given in Fig. (l)(b). Glearly, as k cir- 
cles around (0, 0) or (tt, tt) the direction of b winds twice 
around the unit circle. This "double- winding" explains 
the fact that the degeneracy of the band dispersion is 
lifted "quadratically" as k deviates from (0,0) or (7r,7r). 
On the contrary, the bottom of electron pockets around 
(0,7r) and (tt, 0) are non-degenerate, and b exhibits no 
winding around them. Now let us consider switching on 
a SDW order parameter to nest, say, the fermi surfaces 
around (0, 0) and (tt, 0). Let q be the momentum around 
the (0,0)-fermi surface, and [^^(q)) and |'0(q + (7r,0))) 
be, respectively, the band eigenstates associated with 
the two fermi surfaces. The following matrix element 
A(q) = (V^(q)|M(q)|V^(q+ (7r,0))) determines the SDW 
gap. Here M(q) is a 2 x 2 matrix acting in the or- 
bital space (here we have assumed that after choosing 
a spin quantization axis, the ordered moments lies, say, 
in the =b x-direction). From Fig. (l)(b) one can see that 
H{q) =H{-q) and i7((7r, 0) + q) = i^((7r, 0) - q). How- 
ever, due to the double- winding behavior of b around 
(0, 0) there is a non-trivial Berry phase after k made a 
half circle around the origin, i.e., |?/^(— q)) = —\ip{q)). 
On the other hand the no- winding of b around (tt, 0) im- 
plies |'^(— q+ (7r,0))) = |'^(q+ (7r,0))). Gonsequently if 
M(q) is inversion symmetric, i.e., M(q) = M(— q), we 
have A(q) = — A(— q). Under the assumption that the 
magnetically ordered phase preserves time reversal plus 
a spatial translation (hence A(q) is real), this implies the 
gap function must change sign twice as q moves around 
the fermi surface. Hence there must be two diametri- 
cally opposite nodes. Explicit mean- field calculation [23] 
using the bare H shows that nodes are situated at the 
intersection between the ordering wave vector and the 
(0, 0)-fermi surface. This agrees with the form factor of 
Fig. (3)(bl). 

Next, we come to superconducting pairing. As shown 
in Fig. (3) (a) even for half- filling there are growing in- 
teraction that drives superconductivity! However, these 
interaction are sub-dominant compared with the inter- 
action that promotes antiferromagnetism. Our result 
shows that the two most favorable pairing symmetry 
are an unconventional s (us) and dj;2_y2 like. The u- 
s-wave pairing has s symmetry around the hole fermi 
surface but exhibits s + dj.2_y2 symmetry around the 
electron fermi surface (where the 90 degree rotation sym- 
metry is broken) [8, 20]. It turns out that depending 
on the value of t'^ it is possible for the gap function 
to have nodes on the electron fermi surface. We shall 
return to this point shortly. The form factors fsc.ai}^) 
(A^ca = /5C,a(k)[V^aT(k)V^ai(-k)- T^j]) of thcsc pair- 
ing symmetry are shown in Fig. (3)(b2,b3). Most signif- 
icantly, from our calculation the pairing mechanism can 
he unambiguously determined - the pairing are all driven 
by the inter-pocket pair tunneling[24\ . An example of the 
strongest such process is shown in Fig. (2)(b). As indi- 
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cated by the red and blue arrows, a ±q pair on one fermi 
surface are scattered to a ±k pair on another. Interest- 
ingly, the same scattering also drives the antiferromag- 
netism. This can be seen by noting that the momentum 
transfer between the incoming spin up and outgoing spin 
down electrons is exactly the nesting wavevector (tt, 0) or 
(0, tt) as indicated by the dashed line. Thus the same in- 
teraction also drives antiferromagnetism! One might ask 
would't pairing and SDW interaction require opposite 
sign? No, for the inter band pairing interaction either 
sign will do. This is because the pairing order param- 
eters can choose opposite signs on the two bands thus 
benefit from the positive interaction. At half-filling when 
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FIG. 3: (a)The N=20 RG evolution of the scattering ampli- 
tudes in the SDW and two different types of superconducting 
pairing channels as a function of aln(Ao/A). The parameter 
t[ is set to zero, (bl) The SDW form factor is plotted as the 
momentum moves around the (0,0) or (7r,7r) fermi surface. 
The vertical dash lines are high symmetry directions. The 
curves are obtained by interpolating the data points. Here 
the flat (zero) form factor are associated with the fermi sur- 
faces that are not nested by the ordering wavevector (tt, 0) or 
(0, tt). (b2,b3) The form factor of the two most prominent su- 
perconducting pairing. The most favorable pairing symmetry 
is u-s-wave and the next one is (i^2_„2-wave. 



the nesting is sufficiently good (e.g., when t[ is absent), 
antiferromagnetism overwhelms the superconducting in- 
stability. We propose the reason weak superconductivity 
observed in the stoichiometric compound LaFePO is be- 
cause As^P replacement damages nesting hence allow 
superconductivity to prevail in the competition with an- 
tiferromagnetism. 

The superconducting pairing away from half filling In 
panel (a) and (b) of Fig. (4) we show the N=20 RG evolu- 
tion of the two most favorable singlet and the best triplet 
superconducting scattering amplitudes for x = 0.13. 
The t[ parameter we used here is 0.12 eV, and the bare 
interaction parameters are Ui = 4.0, [/2 = 2.5, J h = 0.7 
eV. Due to the removal of nesting, the antiferromagnetic 
scattering (not shown) is no longer dominant. Like 
half-filling, the most favorable pairing symmetry is the 
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FIG. 4: (a) The N=20 RG evolution of the scattering am- 
plitude associated with the two most favorable singlet chan- 
nel (u-s and (1^2 _y2) and the top triplet channel, (b) The 
form factors of u-s. Note that it changes sign on the electron 
pocket. 



singlet M-s. The pairing mechanism is the inter-pocket 
pair hopping shown in Fig. (2)(b). At t'^ = 0.12 eV the 
form factor changes sign as shown in Fig. (4)(b). As a 
result the superconducting gap has nodes on the electron 
pocket. For all parameters we have studied, the triplet 
pairing channel is never favored. In Fig. (4) (a) we show 
the RG evolution of the best triplet pairing amplitude, 
and it never becomes competitive with the u-s channel. 
We emphasize that while the u-s gap function always 
show a full gap on the hole pockets, it can be gapped 
or gapless on the electron pocket depending on the 
value of t[. For the parameter range we have studied a 
systematic trend is clearly visible: larger t[ makes u-s 
gapless. Under the assumption that our two band model 
describes the band structure adequately, we propose 
that this is the superconducting pairing that has been 
observed experimentally. 

Acknowledgement: We thank Henry Fu for most help- 
ful discussions. DHL was supported by DOE grant num- 
ber DE-AC02-05CH11231. AV was supported by LBNL 
DOE-504108. 



[1] R. Shankar, Rev. Mod. Phys. 66, 129 (1994); J. Polchin- 
ski, Proceedings of 1992 Theoretical Advanced Studies 
Institute in Elementary Particle Physics, edited by J. 
Harvey and J. Polchinski, World Scientific, Singapore 
1993. 

[2] C. Honerkamp, etal Phys. Rev. B. 63, 035109 (2001) 
[3] Y. Kamihara, JAGS, 128, 10012 (2006) and Y. Kami- 

hara, JAGS, 130, 3296 (2008) 
[4] G. F. Ghen, etal arxiv:0803.3790; X. H. Ghen, etal. 

arXiv: 0803.3603; P. Gheng, et.a/ arXiv: 0804.0835 and 

Hai-Hu Wen etal. Europhys. Lett. 82 17009 (2008) 
[5] J. Yang etal. arXiv: 0804.3727 and Zhi-An Ren etal., 

arXiv: 0804.2053 
[6] J. Dong etal. arXiv: 0803.3426, R. H. Liu, etal. 



5 



0804.2105 and Y. Qiu etal. arXiv: 0805.1062; 
[7] C. Cao, P. J. Hirschfeld, H. P. Cheng, arXiv:0803.3236; 
F. Ma, Z. Y. Lu arXiv:0803.3286, D.J. Singh, M.H. Du 
arXiv:0803.0429; K. Kuroki, etal arXiv:0803.3325 
[8] 1. 1. Mazin, etal arXiv: 0803.2740 
[9] V. Cvetkovic and Z. Tesanovic, arXiv: 0804.4678 
[10] H. W. Ou et.al, arXiv: 0803.4328 

[11] C. Cmzar, etal. arXiv:0804.0795; H.-H. Klauss etal 

arXiv:0805.0264 S. Kitao etal arXiv:0805.0041 
[12] K. Ahilan, etal, arXiv:0804.4026 

[13] C. Ren etal. arXiv:0804.1726; F. Hunte etal. 

arXiv:0804.0485 
[14] Y. Nakai, etal. arXiv: 0804.4765 
[15] H. Luetkens etal. arXiv:0804.3115 
[16] L. Shan etal. arXiv:0803.2405 
[17] K. Haule, J. H. Shim, G. Kothar arXiv: 0803.1279 
[18] T. Li, arXiv: 0804.0536; S. Raghu etal. arXiv: 



0804.1113; Q. Han, Y. Chen, Z. D. Wang, Europhys. 

Lett. 82 37007 (2008) 
[19] X. Dai etal. arXiv: 0803.3982, P. A. Lee, X. G. 

Wen, arXiv:0804.1739; X. L. Qi etal. arXiv:0804.4332; 

Z.J. Yao, J. X. Li, Z. D. Wang arXiv:0804.4166; Q. 

Si, E. Abrahams arXiv:0804.2480; G. Baskaran, arXiv: 

0804.1341; Z. Y. Weng, arXiv: 0804.3228; F. J. Ma, 

Z.-Y Lu, T. Xiang, arXiv:0804.3370 and K. Seo, B. A. 

Bernevig and J. P. Hu, arXiv:0805.2958. 
[20] Z. H. Wang etal. arXiv: 0805.0736, Y. Wan, Q. H. Wang 

arXiv: 0805.0923 
[21] Fa Wang et al, arXiv:0807.0498. 
[22] M. Berry, Proc. R. Soc. London, A 392, 45 (1984) 
[23] Ying Ran et al, arXiv:0805.3535. 

[24] H. Suhl, B.T. Matthias, L.R. Walker, Phys. Rev. Lett. 
12, 552 (1959). 



